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The compression of the layered carbon nitride C 6 N 9 H 3 -HC1 was studied experimentally and with density 
functional theory (DFT) methods. This material has a polytriazine imide structure with CI" ions contained 
within C 12 Ni2 voids in the layers. The data indicate the onset of layer buckling accompanied by movement of 
the CI" ions out of the planes beginning above 10-20 GPa followed by an abrupt change in the diffraction 
pattern and c axis spacing associated with formation of a new interlayer bonded phase. The transition 
pressure is calculated to be 47 GPa for the ideal structures. The new material has mixed sp 2 -sp 3 
hybridization among the C and N atoms and it provides the first example of a pillared-layered carbon nitride 
material that combines the functional properties of the graphitic-like form with improved mechanical 
strength. Similar behavior is predicted to occur for Cl-free structures at lower pressures. 
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High pressure investigations of carbon nitride materials have been stimulated by theoretical predictions that 
dense sp 3 -bonded C 3 N 4 phases would display low compressibility and high hardness values 1,2 . Despite 
many synthesis attempts including use of high pressure the formation of crystalline C 3 N 4 polymorphs 
remains an elusive goal. However an sp 3 -bonded carbon nitride imide C 2 N 3 H with a defective wurtzite structure 
has been produced in laser-heated diamond anvil cell (DAC) experiments and this material was recoverable to 
ambient conditions 3 ' 4 . Low density sp 2 bonded polymeric or layered graphitic carbon nitride materials (g-CNMs) 
are also well known 5 and these have attracted attention as metal-free redox catalysts, photocatalysts and electro- 
ceramic materials 6 ' 7 . They also provide precursors for investigation of pressure-induced transformations into 
dense sp 3 bonded phases. We investigated the room temperature compression behavior of one such well- 
characterized g-CNM of composition C 6 N 9 H 3 *HC1 8 to 70 GPa using synchrotron X-ray diffraction experiments 
in a diamond anvil cell (DAC) combined with density functional theory (DFT) calculations carried out to 
100 GPa. The combined experimental and theoretical results indicate the initial onset of layer buckling and 
movement of the CI" ions out of their sites within the planes followed by a phase transformation into a structure 
containing interlayer C-N bonds between sp 3 hybridized atoms. The new material constitutes a new example of a 
pillared-layered gCNM with mixed sp 2 -sp 3 bonding. Related bonding changes have been recorded in graphitic C 
and BN materials at high pressure 9 " 14 . 

A large class of oligomeric, polymeric and graphitic carbon nitride materials are determined to have structures 
based on heptazine (tri-s-triazine, C 6 N 7 ) rings linked by -N= or -NH- groups 15 . A second structural motif based 
on polymerized triazine (C 3 N 3 ) units is found among graphitic nanocrystalline C 3 N 4 prepared by chemical vapor 
deposition 16 . Triazine units are also the building blocks in bulk crystalline compounds synthesized from pre- 
cursors including melamine (C 3 N 6 H 6 ) and cyanuric chloride (C 3 N 3 C1 3 ) under high-P,T conditions (0.5-3 GPa, 
300-500°C). Further studies have obtained crystalline graphitic polytriazine imides (PTI) with related structures 
by reactions in molten salt mixtures (e.g. KCl-LiCl eutectics) 517 . The s-triazine rings linked by bridging -N= or - 
NH- units lead to large voids within the graphitic layers. In the case of C 6 N 9 H 3 *HC1 studied here the Cl~ ions 
occupy the centers of the voids at ambient pressure. The layers are stacked to provide an alternating pattern of 
triazine and Ci 2 Ni 2 rings along the c axis (Fig. 1). 

Results 

Experimental findings. The starting compound g-C 6 N 9 H 3 *HCl is descibed by a hexagonal unit cell with P6 3 /m 
symmetry (Fig. I) 818 . The X-ray patterns at low pressure (Fig. 2) are dominated by the (002) interlayer reflection at 
20 ~ 8° that is indicative of a layered graphite-like structure with an ABAB stacking along the c axis. The 
polytriazine imide layers contain C 12 N 12 H 3 voids hosting CI" ions derived from the synthesis reaction 
between melamine and cyanuric chloride, and additional H + ions are bound to the N atoms of the triazine 
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Figure 1 | Structural features of graphitic C 6 N 9 H3-HC1 at ambient pressure, (left) Top view of one plane of the structure showing the triazine (C 3 N 3 ) 
rings connected by -NH- groups to form large C 12 N 12 voids that are occupied by CI" ions. The accompanying extra H + ion is attached to one of six 
possible N positions from the triazine units surrounding the large ring 616 . The H atoms have been omitted for clarity, (right) View down c of two 
adjacent layers of the structure showing the ABAB stacking sequence that places one triazine ring above the C 12 N 12 void in successive planes. The H and CI 
species are not shown for clarity. 



rings 818 . Only one of six available N sites is protonated in this way so 
that the P6 3 /m space group used to analyze the data represents a 
spatially averaged solution. The ABAB stacking of the graphitic 
layers creates a structure in which one half of the triazine rings 
within each layer is positioned above or below C 12 N 12 voids in 
adjacent layers, while the others overlap triazine rings in the layers 
above and below (Fig. 1). 

The prominent (002) interlayer reflection shifts rapidly to smaller 
d spacing (20 ~ 10°) with increasing pressure (Fig. 2). Above 10 GPa 
the diffraction peaks become significantly broadened but the general 
features of the graphitic material are still recognizable up to P ~ 
36 GPa. However at P = 40 GPa the nature of the pattern has chan- 
ged to become dominated by a main broad asymmetric peak near 13° 
20. This result signals a structural change into a new high-density 
form (Fig. 2). After the transformation has occurred the peak posi- 
tions vary little with continued compression up to 70 GPa indicating 
that the high pressure structure is significantly less compressible than 
the graphitic layered phase. The substantial peak broadening that 
occurred as a function of pressure did not permit us to carry out 
detailed refinement of the diffraction data. Instead we investigated 
the nature of the structural changes at high pressure using density 
functional theory (DFT) calculations, that enabled us to gain under- 
standing of the pressure induced transformation by comparing the 
predicted diffraction profiles with experiment (Fig. 3). 

DFT calculations. In a first series of DFT studies we geometry- 
optimized the structure of g-C 6 N 9 H3*HCl between P = 0- 
100 GPa in 10 GPa steps, starting each optimization from the 
ambient pressure lattice spacings obtained from experiment (a,b = 
8.42 A; c = 6.34 A). We generated a proton ordered version of the 40- 
atom crystallographic unit cell within space group PI, removing all 
local symmetry constraints to allow the atoms to relax both within 
and away from the graphitic planes. This differs from our experi- 
mental data analysis that assumed a spatially averaged P6 3 /m unit 
cell 818 . The proton-ordered structure exhibits a small shear distortion 
between adjacent planes along (010) to maximize electrostatic 
interactions between the H + and Cl~ ions within layer voids. This 
leads to a small separation between a and b cell parameters that 
increases slightly with compression but that remains insignificant 
until above 20 GPa (Fig. 4). At P = 0 GPa the optimized lattice 



parameters were found to be a,b = 8.37 A; c = 6.23 A, slightly 
smaller than those found experimentally (8.42, 6.34A) as is typical 
for LDA results. During compression there is a rapid decrease in the 
c-axis parameter consistent with experiment: between 0-40 GPa 
there is a 17% reduction in c, compared with 3.7% and 5.8% in the 
a and b lattice parameters (Fig. 4). Between 60-70 GPa a sudden 
collapse of the c-axis occurred in these initial compression experi- 
ments, and then this parameter remained approximately constant up 
to 100 GPa (Fig. 4). 

These initial results showed that a dramatic change had occurred 
in the structure as the c axis dimension was reduced by compression, 
implying a separation between the initially planar layers of approxi- 
mately 2.0 A, without considering the reduction in distances due to 
layer buckling. Such a small interlayer separation could lead to the 
onset of C-N interlayer bonding and we investigated this possibility 
by analyzing the fully relaxed structures as a function of densifica- 
tion. The results led us to identify the structure optimised at high 
pressure as a new phase of the material, characterized by the occur- 
rence of bonding between a subset of sp 3 hybridized C and N atoms in 
adjacent layers (Fig. 5). This new phase is termed the ILB1 structure 
to distinguish it from the densified graphite-like phase with buckled 
but independent layers obtained at lower compression (ILB0). The 
C-N bonding occurs between atoms associated with triazine rings 
that in the ambient pressure phase overlapped each other in adjacent 
layers; the others located above and below the C 12 N 12 voids remain 
sp 2 -hybridised giving rise to a unique new pillared-layered structure 
that contains a mixture of sp 2 - and sp 3 -bonded C and N atoms 
(Fig. 5). 

Once we had identified the occurrence and connectivity of the new 
high-pressure phase, we examined the relative enthalpies of the ILB0 
and ILB1 polymorphs over a wide pressure range, starting the studies 
using the compressed but still graphite-like (ILB0) structure obtained 
at 60 GPa and then gradually increasing P, compared with results 
during gradual decompression of the ILB1 phase obtained at 70 GPa 
(Fig. 6). Each calculation was initiated from the geometry- optimized 
structure at the immediately adjacent pressure value. The C-N con- 
nectivity was maintained throughout in both series of calculations, 
indicating that the ILB0 and ILB1 phases correspond to separate 
functions in the potential energy surface of the system and that their 
interconversion is an activated process that does involve interlayer 
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Figure 2 | Angle dispersive synchrotron X-ray diffraction data for 
C 6 N 9 H 3 -HC1 obtained up to P = 70 GPa (A 0 = 0.4441 A). Peaks of 
interest are labelled with hkl designations according to the P 63/01 space 
group, (a) Selected patterns showing the onset of peak broadening for the 
layered phase above 10 GPa followed by the rapid emergence of the new 
pattern above 40 GPa. (b) Variation of d spacings measured for the 
principal peak maxima assigned to the low pressure graphitic phase as a 
function of pressure. 

but not intralayer C-N bond breaking or formation. Comparing the 
calculated enthalpies we identified 47 GPa as the equilibrium trans- 
ition pressure between the ILB0 and ILB1 polymorphs at T = 0 K 
(Fig. 6). Following the compression studies we relaxed the ILB1 
structure to ambient pressure, and found that the pillared connec- 
tivity was retained. Similar calculations repeated for the hypothetical 
Cl-free compound showed a similar structural transformation occur- 
ring at high pressure, but with the isoenthalpic point between ILB0 
and ILB1 phases reduced to ~5 GPa. That result indicates that vary- 
ing the CI content of these gCNM phases should cause the transition 
pressures and recoverability of the pillared-layered material to be 
adjusted experimentally. 

In order to fully characterize the structural changes associated 
with compression of C 6 N 9 H 3 • HC1 from its ambient pressure graph- 
itic phase to the denser ILB1 polymorph at high pressure, we mea- 
sured a set of dihedral angles in the minimum energy structures 
relaxed at each pressure. These angles (Fig. 7a) are chosen to describe 
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Figure 3 | Comparison of experimental (black line) with the DFT 
predicted (ILB0 phase, red line) X-ray patterns for the graphitic phase at 
various pressures. At 40 GPa the blue pattern corresponds to dense forms 
with interlayer bonding (ILB1) as described in the text. 
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Figure 4 | Changes in the DFT calculated lattice parameters for a 
proton ordered version of the g-QNgPVHCl structure as a function of 
pressure. 
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Figure 5 | Left: The ILB1 phase of C 6 N 9 H 3 , HC1 obtained at high pressure from DFT calculations. In the real structure studied experimentally the 
Cl~ displacements above and below the layers can occur in a random or locally ordered fashion. Right: picture highlighting the formation of pillars 
with three-dimensional connectivity between the carbon nitride layers produced by interlayer bonding within a mixed sp 2 -sp 3 bonded state; the pillars are 
connected within the planes by relatively undistorted triazine rings. 



the distortion of the initially planar graphitic layers: 1) the internal 
dihedral angle of a protonated C 3 N 3 ring (ABCD); 2) the dihedral 
angle formed by the skeletal proton (i.e. the bridging NH group) with 
respect to the neighbouring triazine ring (ABFG); and 3) the internal 
dihedral angle of a non-protonated C 3 N 3 ring (HIJK). All these 
values are identically zero for fully planar layers and any deviations 
thus record buckling of the gCNM framework. Achieving a value of 
60° in any of these parameters indicates a change of hybridization at 
the central C and N atoms from sp 2 to sp 3 . At each pressure reported 
in Fig. 7 we employ the thermodynamically stable phase from the 
DFT calculations, i.e. ILBO below the transition pressure of 47 GPa, 
and ILB1 above. 

At ambient pressure the ILBO structure is nearly planar with all 
dihedral angles <2°, as expected for a graphitic phase with deloca- 
lized sp 2 bonds. However, this condition is lost as soon as the struc- 
ture is compressed even slightly. By 10 GPa the dihedral angles HIJK 
and ABFG differ noticeably from zero but those within the C 3 N 3 
rings (ABCD) remain constant. This reveals a buckling of the layers 
around the bridging -NH- groups (ABFG) beginning at low pressure 
that is accompanied by a displacement of the Cl~ ions out of the 
graphitic planes (Fig. 7b). Between 40-50 GPa the ABCD dihedral 
angle jumps from 2° to 65°, indicating the sp 2 to sp 3 hybridization 
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Figure 6 | Calculated enthalpy as a function of P for the layered-pillared 
ILB1 phase relative to the layered graphitic (ILBO) polymorph of 
C 6 N 9 H 3 • HCI and of the hypothetical Cl-free C 6 N 9 H 3 system. 
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structural phase transformation occurs at 47 GPa. (b) Displacement of the 
Cl~ ions along the c axis relative to the initial graphitic planes within the 
calculated ILBO structure. 
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change that has taken place around the C and N atoms associated 
with the triazine rings. At this point the reduction in the c axis is 
such that the interlayer C-N bonding has occurred to yield the ILB1 
structure. The structural changes can also be correlated with the 
behavior of the CI" ions. While these remain within the layers, they 
are clamped tightly by strong H-bonding interactions with the NH 
groups around the Ci 2 Ni 2 ring sites. On application of pressure, 
these voids are diminished in size and they become too small to host 
the CI" ions that are progressively squeezed out of the plane. This 
effect becomes more marked with increasing densification until it 
reaches nearly 0.6 A by 70 GPa. These CI" displacements also lead to 
structural disorder in the real system, as the CI" ions can be located 
either above or below the planes and this effect can contribute to the 
peak broadening observed in the experimental X-ray diffraction pat- 
terns. Once removed from the graphitic layers and into the interlayer 
space the CI" ions also become free to be displaced along the (a,b) 
directions. 

Comparison with experiment: ILBO vs ILB1 formation at high 
pressure. Following our theoretical prediction of formation of a 
novel ILB1 structured phase at high pressure, we compared our 
calculated X-ray diffraction patterns for this material with experi- 
mental data (Fig. 8). The diffraction pattern recorded experimentally 
at 70 GPa showed broad features extending throughout the 5-18° 26 
range, that could not be refined to provide structural details. 
However comparison with the calculated results for the predicted 
ILB1 and ILBO structures at the same pressure is revealing. The 
predicted X-ray pattern for the optimized ILB1 phase at 70 GPa 
matches the observed diffraction profile, with the main diffraction 
intensity concentrated near 13° 26, along with other features 
occurring near 10° and 6.5° 26. The pattern is significantly diffe- 
rent from that expected for the ILBO structure at the same 
pressure, leading us to conclude that the high pressure structural 
transformation observed experimentally corresponds to a transi- 
tion between ILBO and ILB1 phases. 

Low pressure V(P) results. The computational results enabled us 
not only to identify the structural features of the new high-pressure 
ILB1 polymorph of the carbon nitride material under study, but also 
to further characterize the evolution of the graphitic ILBO phase at 
intermediate pressures, highlighting the occurrence of buckling 
distortions and displacements of the CI" ions. We carried out a 
detailed analysis of our angle dispersive X-ray diffraction data 
below 36 GPa to study the compression behavior of the ILBO 
structure and validate the computational results in relation to the 
experimental measurements. Refinements were performed within 
P6 3 /m using PowderCell® software by first assigning the position 
of the main diffraction peak to the (002) reflection of the ABAB 
layer- ordered structure, using the DFT results to guide the 
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Figure 8 | Comparison of the experimental XRD pattern obtained at 
70 GPa (top) vs simulated patterns generated from DFT calculations 
(bottom), derived from geometry optimised ILB1 and ILBO structures of 
C 6 N 9 H 3 • HC1, with or without interlayer bonding occurring at the same 
pressure. 



structural analysis in the highest pressure ranges where substantial 
peak broadening was present (Fig. 2). The variations in refined a,b,c 
cell parameters led to the V(P) plot shown in Figure 9. 

Between 20-30 GPa we note a slight deviation of the V(P) relation 
upwards from the expected equation of state (EoS) curve obtained by 
analyzing the data using the spatially averaged space group P6 3 /m. 
This can occur due to layer buckling along with displacements of the 
CI" ions away from the graphitic planes, that become significant 
above 20 GPa as shown by the DFT results (Fig. 7b). We attempted 
to analyze our data in this range using lower symmetry space group 
solutions P6 3 , P6 and PI but these all gave similar R wp and % 2 values 
despite the increased structural degrees of freedom. 

We fitted the V(P) data up to 21 GPa using a third order Birch - 
Murnaghan (BM-3) equation of state: 

P( V) = 3K 0 f( 1 + 2ff'4l + I (K' 0 - 4)/ 



with K 0 the zero pressure bulk modulus and K 0 its pressure deriv- 
ative, and the volume strain variable expressed by 



1 



This fit yielded K 0 = 32.8 ± 7.3 GPa and K 0 ' = 4.8 ± 0.5, in good 
agreement with our DFT results for the ILBO polymorph (K 0 = 
29.1 GPa; K 0 ' = 6.7). The compressibility is comparable with crys- 
talline graphite (K 0 = 34 GPa), although poly crystalline samples 
have been reported with K 0 values extending down 10 GPa 19 . 

Fitting our DFT V(P) results for the high-pressure ILB1 
QNgL^'HCl structure to a BM-3 EoS yields a bulk modulus K 0 = 
130 ± 5 GPa and K 0 ' = 3.57 ± 0.2, clearly showing that the ILB1 
structure has much lower compressibility than the sp 2 bonded graph- 
ite-like ILBO phase (K 0 = 29.1 ± 4.1 GPa). The new high pressure 
ILB1 phase, despite partly retaining the void intralayer structure 
of the low pressure ILBO polymorph, has therefore substantially 
higher structural strength, a feature of potential interest for practical 
applications. 

Discussion 

We conducted a high pressure study of graphitic C 6 N 9 H3'HC1 using 
complementary experimental and theoretical approaches. The 
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Figure 9 | V(P) plot obtained from our experimental data (circles) for 
g-C 6 N 9 H 3 .HCl analyzed within the spatially averaged space group P6 3 /m. 
The data were fitted up to 2 1 GPa using a third order Birch-Murnaghan 
equation of state (K 0 = 32.8 ± 7.3 GPa and K 0 ' = 4.8 ± 0.5) and 
extrapolated to 36 GPa (line). 
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graphitic layered compound C 6 N 9 H 3 • HC1 was subjected to pressures 
of up to 70 GPa in a diamond anvil cell, and its structural behaviour 
was examined using synchrotron X-ray diffraction. The structural 
model proposed by Zhang et al. with a layered hexagonal (P6 3 /m) 
structure was employed for the experimental data analysis, for which 
it provides a useful description up to a pressure of —40 GPa. Close 
examination of the V(P) relation indicates an anomaly near 20 GPa 
that can be explained by a deviation of the structure from P6 3 /m 
symmetry that can be interpreted from the computational results. On 
increasing pressure from ambient conditions, we observe three dis- 
tinct structural rearrangements: at ambient pressure the carbon 
nitride layers are atomically flat and host the CI" ions within the 
layers; between 10-20 GPa the layers buckle and the CI" ions dis- 
place vertically but not horizontally with respect to the layers; above 
20 GPa the CI" ions can displace significantly both horizontally and 
vertically, making the P6 3 /m model less and less appropriate for a 
detailed structural description. 

The calculations identify the formation of a new pillared structure 
(referred to as ILB1), that is thermodynamically stable compared to 
the graphitic phase above 47 GPa. Simulated XRD patterns using 
graphitic and pillared phases from the computational work were 
compared with the experimental XRD data, and these indicate that 
the predicted pillared structure provides the most plausible match to 
the observed data above 40 GPa. This pillared structure is formed by 
buckled sp 2 -hybridized platelets joined by sp 3 -hybridised "tubes" 
running along the (001) direction (Fig. 5). The close match between 
the experimental data and that of the ILB1 structure suggests that 
the initially graphite-like system undergoes a change in hybridi- 
zation and experiences interlayer C-N bonding as predicted by the 
calculations. 

In conclusion, our experimental and DFT results indicate that 
compression of the layered compound g-C 6 N 9 H 3 *HCl at ambient 
temperature results in a polymorphic transformation leading to a 
new interlayer bonded phase. This has mixed sp 2 -sp 3 bonding 
between C and N atoms in adjacent layers and it gives rise to a 
new family of pillared-layered compounds that can combine the 
essential characteristics of the 2D material with buckled C x N y (H z ) 
planes and the mechanical properties of a 3D-bonded compound. 

Methods 

The starting compound g-C 6 N 9 H 3 *HCl was formed from a 1 : 2 mixture of melamine 
(C 3 N 6 H 6 ) and cyanuric chloride (C 3 N 3 C1 3 ) reacted under high-P,T conditions 
(1.5 GPa, 550°C) in a piston-cylinder device. The synthesis details and structural 
characterization of this material were described previously 818 . 

Synchrotron X-ray diffraction studies were carried out at up to 70 GPa using angle 
dispersive techniques at ESRF ID27 (Grenoble, France: l Q = 0.3738 A: Figure 2). 
Additional data were obtained at Diamond Light Source 115, UK (1 Q = 0.4441 A). 
Samples were loaded into mechanically (four post Mao-type symmetrical cells) or gas- 
driven He membrane diamond anvil cells (DACs). Diamonds with 200-300 um culet 
sizes were used to collect data up to 70 GPa. Re gaskets pre-indented to —30 um with 
gasket holes —80 um in diameter were drilled by electro -erosion. Samples were packed 
into the cell chamber inside a glove box without a pressure-transmitting medium 
(PTM). Ruby fluorescence spectra indicated that the pressurisation environment 
remained reasonably "hydrostatic" throughout most of the runs: pressure gradients of 
only —1-2 GPa across the cell were detected and the ruby signals showed little 
broadening. The samples were loaded without a PTM to avoid interference with the 
weak and often broadened diffraction signals from the low-Z carbon nitride material, 
and to avoid potential effects associated with intercalation of fluid species such as He, N 2 
or H 2 into the layered structure. X-ray diffraction data were recorded using MAR image 
plate or CCD detectors and analysed using Fit2D 20 , FullProf 2 1 and GSAS 22 ' 23 software. 

Electronic structure calculations were carried out under periodic boundary con- 
ditions within the local density approximation (LDA) of density functional theory 
(DFT) using the planewave code CASTEP 24 . We used ultra-soft pseudopotentials and 
chose the cut-off energy (450 eV) and k-point grid size (3 X 3 X 4) to obtain best 
converged results 25 . Calculations were performed under hydrostatic pressure condi- 
tions: full geometry optimizations were carried out at each pressure within space 
group PI using a 40-atom unit crystallographic cell. Prior to each geometry 
optimization, the atoms were subjected to small random displacements to produce 
symmetry-free initial structures. 
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